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Abstract 

We investigate a propagation of solitons for nonlinear Schrodinger 
equation under small driving force. The driving force passes through 
the resonance. The process of scattering on the resonance leads to chang- 
ing of number of solitons. After the resonance the number of solitons 
depends on the amplitude of the driving force. The analytical results 
were obtained by WKB and matching method. We bring two examples 
of numeric simulations for verifying obtained analytical formulas. 

Introduction 

Nonlinear Schrodinger equation (NLSE) is a mathematical model for wide 
class of wave phenomenons from the signal propagation into optical fibre ^ EJ 
to the surface wave propagation [3] . This equation is integrable by inverse scat- 
tering transform method 4_ and can be considered as an ideal model equation. 
The perturbations of this ideal model lead to nonintegrable equations. Here 
we consider such nonintegrable example which is NLSE perturbed by the small 
driving force. 

ia t * + 9^+|*| 2 * = e 2 /e iS/£2 , 0<£<1. (1) 

The perturbed NLSE in form JIJ is relative to describing of nonlinear effects 
of optical soliton propagation in the presence of an input fast oscillating forcing 
beam [5]. 

The most known class of the solutions of NLSE is solitons 0] . The structure 
of this kind of solutions is not changed in a case of nonperturbed NLSE. The 
perturbations usually lead to modulation of parameters of solitons [El Q, for 
driven NLSE see also |Hj- Number of solitons does not change. 

In this work we investigate a new effect called the scattering of solitons on the 
local resonance. Typical picture for this process may be obtained by numerical 
simulations (see fig. 1) . 
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Let us explain the terms 'scattering' and 'resonance' which is used in the 
work. Usually one says 'scattering' for process when a wave is changed after ob- 
stacle. The same process is studied in our work for solitons which are nonlinear 
waves of NLSE. The obstacle is a line where the external force has a resonance in 
the equation for first-order term in the asymptotic solution for perturbed NLSE. 
Under 'resonance' we understand standard phenomenon of grows solution be- 
cause of oscillating external force. The resonance phenomenon takes place in a 
thin domain near some curve usually called a resonance curve. So we use the 
term 'local resonance'. 

We consider the process of scattering in detail and obtain the connection 
formula between pre-resonance and post-resonance solutions. In general case 
the passage through the resonance leads to changing of the number of solitons. 
This effect is based on the phenomenon of soliton generation due to the passage 
through the resonance by the external driving force [H]- 

We found that the scattering of solitary waves on the resonance is a general 
effect for the wave propagation in a nonlinear media. In this work we investigate 
this effect for the simplest model. It allows to show the essence of this effect 
without unnecessary details. 

This paper has the following structure. The first section contains the state- 
ment of the problem and the main result. The sections 2-4 are reviewed asymp- 
totic analysis of the problem which was developed in ^HI- The second section 
contains the asymptotic construction in the pre-resonance domain. In the third 
section we construct the asymptotic solution in the neighborhood of the res- 
onance curve. The fourth section of the paper is devoted to construction of 
the post-resonance asymptotics. Asymptotics are constructed by multiple scale 
method an( i matched At last five section contains the results of nu- 

merical simulation which justify the obtained asymptotic formulas. 

1 Statement of the problem and result 

We study the perturbed NLS equation JJJ with a special phase of the driving 
force: S/e 2 = e 2 t 2 . The amplitude / = f(ex) is a smooth and rapidly vanished 
function. The small parameter in the right hand side of the equation is defined 
as e 2 only for a convenience. 

The goal of our work is the study of a slow evolution of the solution with a 
small amplitude for the equation QJ. In the general case the small amplitude 
solutions are defined by the linear Schrodinger equation. However there exists 
a magic relation between scales of independent variables of a carrier wave and 
an amplitude of its envelope function which leads to the nonlinear Schrodinger 
equation for the envelope function of the leading-order term of the asymptotics. 
This relation was observed for different physical problems in pioneer works 1 
El E| • An example of such relation give the following substitution 

^ = eU{f-2,xi,e), t 2 = e 2 t 1 X\=ex, 
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then the equation Q has the form: 

id tl U + d* 1 U+\U\ 2 U = s- 1 fe iS / s2 . (2) 

The strong perturbation with rapid phase for NLSE may be considered as 
a model for the high-frequency heating of plasma an d leads to the phe- 
nomenon of the scattering of solitons. 

It is known if / = then there exists a soliton solution of the NLSE. 
We'll show that for / ^ the number of solitons may be changed due to 
passage through the resonant curve. The resonant curve defines as a line where 
the frequency of the driving force is equal to the cigcnfrcquency of linearized 
Schrodinger equation. In our case this curve is t = 0. After passage the reso- 
nance the number of solitons depends on amplitude of the perturbation on the 
resonance curve t = 0. 

In the simplest case the phase is linear function with respect to t S/e 2 — 
ut, to — const. In general situation the constant frequency of the driving force 
does not lead to the scattering of solitons. Let us investigate the driving force 
with slowly varying frequency. The most simplest dependence on t for lj has 
a form to = e 2 t/2. Hence S = e 4 t 2 /2. Namely the equation with this form of 
the phase function of the driving force is studied in this work. The amplitude 
/ of the driving force admit an additional dependency on e 2 t but it leads to 
complicated formulas and no more. 

Let us formulate the result of this work. Below we use the following variables 
Xj — — x , tj — s,^ t j j — _L . . . 

Let the asymptotic solution of Q be 

^f(x,t,e) = e u (xi,t 2 ) + 0(e 2 ) as - c < t 2 < 0, 
where c = const > and u (xx, t 2 ) satisfies 

d t2 u +d 2 xi u +| u | 2 u= 

and initial condition 

u \t 2 =t° = hi(xi), t° — const < 0. 
Then in the domain < t 2 < c the asymptotic solution of has a form 

*(x,t,e) = e v (x 1 ,t 2 ) + 0(e 2 ), (3) 

where v (xi, t 2 ) is a solution of NLSE with initial condition 

v\ t2 =o=u(xi,0) + (l-i)y/7rf(x l ). (4) 

Formula Q is the main result of the paper. This formula connects the main 
order term u of the asymptotic solution before the scattering, external driving 
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force / and initial condition v for the main order term of the asymptotic solution 
after the scattering. This formula is derived in the end of section 4. 

Let us explain the result for soliton solution. If in the domain — c < t% < 
the solution has iV-soliton form then in the domain < t 2 < c the number of 
solitons is defined by initial condition |4"|. 

Given analysis is valid at \t 2 \ < C, for VC = const. When \e 2 t\ ^> 1 the 
perturbation force will modulate the parameters of the soliton solution, sec 

mm 



2 Incident wave 

In this section we construct the asymptotic solution of equation Q in pre- 
resonance domain. This solution contains two parts. The first part is a specific 
solution of the nonhomogeneous equation. This solution oscillates with the fre- 
quency of the driving force. The amplitudes are determined by an algebraic 
equations. The second part of the solution is a solution of the homogeneous 
equation. The solution contains an undefined function due to integration. This 
undefined function usually determines by initial condition for the Cauchy prob- 
lem. 

We construct the formal asymptotic solution in the WKB-like form 

#(x,t,E) = eu (x 1 ,t 2 ) + e 3 u (xi,t 2 ) + e 2 B {xi,t 2 )exp(iS/e 2 ) + 

£ 4 (^Bi (x u t 2 )exp(iS/s 2 )+ B-i (x u t 2 )exp(~iS/e 2 )^ (5) 

+e 5 B 2 {x u t 2 )exp(2iS/e 2 ). 

To determine the coefficients of the asymptotics substitute (JSJ into equation 
CJ. It yields 



e 2 i - S' B - fj exp{iS/e 2 ) + e 3 \i u t2 + u XlXl +\ u \ 2 u 
+e 4 ^ - S' Bi +i Bt 2 + B XlXl +2| u \ 2 B j cxp{iS/e 2 ) + 

S' B-i + u 2 B *\ cxp{-iS/e 2 )^j 



e b ( 1 u t2 + u XlXl +2| u\u + u u +2u\b\ + 
25' B 2 + u * B exp{2iS/e 2 )^j = e Q R(t 2lXl ). 
The residue part of the asymptotics has a form 

R(t 2 , Xl ) = 0(\ 1 | 3 +e 3 | u | 3 +e 6 | Bi \ 3 + e 6 \ B-i \ 3 + e 9 \ Bi | 3 ). (6) 
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Collect the terms with the same order of e up to the order of e 5 and reduce 

1 3 

similar terms. It yields differential equations for u, u and algebraic equations 

2 4 5 



for B, B±i and B2- 



i u t2 + u XlX1 +| u | 2 u= 0, (7) 



i u t2 + u XlXl +2\ u | 2 u + u 2 u * = -2| B | 2 u, (8) 

- S' B= f, (9) 

S' Bi= i B t2 + B XlXl + \ u | 2 B, (10) 

S' B-i=-u 2 B*, (11) 

~2S'B2=u*B 2 . (12) 

1 3 

The functions w, u are uniquely determined by initial conditions at the mo- 
ment t2 =t°. We suppose that t° = const < and 

1 3 

u | i2=t o = fti(ari); u | t2=t o = ft 3 (a;i); 

where functions ft-i, are smooth and rapidly vanish as |xi| — > ±00. 

1 3 

It's known the solutions U and u exist for bounded values of ti-, see [141 115) . 

3 

Remark The solution of the equation for u contains growing terms as t% — > 
00, see for example JSl- These terms are secular as t<i ~ e -1 . But we do not 
consider such long times in this work. 

The coefficients of the representation © have a singularity as S' — > 0. The 

3 

order of singularity of Bk is easy calculated. 

B=0(«- x ), Bi=0(t~ 3 ). 
3 

To determine the asymptotics of u as £2 —> —0 we construct the solution of 
the form 

u= u (_1 ' 0) (asi,fe) + In |i 2 | u (0,1) (a;i, *a) + *2 hi |t 2 | u (M) (zi, t 2 ) + u{ Xl M)- 

(13) 

Substitute this representation into equation (jHJ and collect the terms of the same 
order with respect to t% . It yields equations for coefficients of the asymptotics 



3 
it 



(-^> = J 2|/| 2 i, 
= -iL(u (- 1 - )), 

L(«) = i* 2 In |i a |£(u (1A) ) + i I (14) 
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Here L(u) is a linear operator of the form 

L(u) = id t2 u + dl x u + 2\u | 2 u+ u 2 u* . 

Functions and u are determined from algebraic equations. 

These functions are bounded as —const < t<i < 0, const > 0. 

3 

The function u is a solution of nonhomogeneous linearized Schrodinger equa- 
tion. The right hand side of the equation is a smooth function as —const < 
^2 < 0, const > 0. The solution of this equation can be obtained using re- 
sults of JS]. In particularly if u is N-solitons solution of NLSE then exists the 
bounded solution of nonhomogeneous linearized Schrodinger equation 1141) as 
—const < ti < 0, const > 0. 

The asymptotic form (J5J allows to solve equation JIJ up to the order £ 6 . 
To obtain more accurate approximation one have to include terms without fast 
oscillating of the order o(e 3 ) into asymptotic solution (JjJJ. Therefore we define 
the domain of validity of (J5J by following relation: 

e 6 R(t 2 , Xl ) = o(e 3 ), e^0. 

Coefficients of (jSJ have singularity at ti = 0. The residue part increases as 
ti — > 0. From formulas 10 and l(7|l- H12|) one can easily obtain the behaviour of 
the residue part: 

R{t 2 , Xl ) = 0{t- 3 +e e t- 9 ), t 2 ^-0. 
It yields the domain of validity of 10 

-t 2 >£ or -t>e _1 . 

3 Scattering 

In the neighborhood of the point t 2 = the frequency of the driving force 
becomes resonant. Formally it means representation @ is not valid. 

In this part of the work we construct another representation for the solution 
of equation . This representation is valid in the neighborhood of the resonance 
line t 2 = 0. 

1 2 

^>(x,t,s) = £ w (xi,h)) +£ 2 w (xi,ti) + 

3 1 3 
£ 3 hl£ w (xi, t\) + £ 3 W (xi, ti) £^0. (15) 

Here we use a new scaled variable t\ — t 2 /e. Representation l|15|l is matched 
with (JSJ. It means these formulas are equivalent up to value o(e 5 ) as t 2 — * —0. 
The coefficients of (|15|) are determined by ordinary differential equations (|16|) . 
ljTH|) . lf^U|) and matching conditions. 
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To obtain the behaviour of the coefficients of (|f 5|) as t\ — > — oo match <|1 511 
with (0. Write © in terms of ti 



e 2 



(x, t,e)=e(l{x 1 ,0)- (t^f + it^f) exp(iS/e 2 )^j + 

(d t2 u {x 1 ,t 2 )\ t2=0 t 1 + t-h\f\ 2 1 (x u o) + o(t- 2 )^j + 

e 3 \ne( - iL(2i\f\ 2 u)| ta=0 + o(l)] + 



+ £3 {\ d l u(x 1 ,t 2 )\t2=otl+u(x u 0) + o{l)\, l<-ii<£ \ e^O. 

To obtain equations for coefficients of l|15|) substitute l|15f) into equation ■ 
It yields 

e 2 {[d tl w -fexp(iS/e 2 )^J + e 3 (d tl w +d 2 xi w + 7 | w \ w \ + 

(3 2 12 1 2 \ 

9 tl u> w++w 2 w*+2 1 \w\w ) = e 5 p(h,xi,e). 

The function p(ti,a;i,e) can be represented in the form 

p(t 1: x 1: e) = 0(| w\ 2 w +d 2 Xl w +e\ w | 3 + e 4 | u> | 3 ). 

Collect the terms of the same order with respect to e. As result we obtain the 
equations for coefficients of (|15|) . 

id tl i=/exp(^/2). (16) 

The matching conditions give (xi,0) £i — > — oo. The solution of this 

problem is represented in terms of Fresnel integral 

i i r' 1 

w=u(x 1 ,0)-if(x 1 ) exp(i6 2 /2)d8. (17) 

J — oo 

Equations for higher-order terms are 

id tl w= -d 2 Xl w+\w\ 2 w, (18) 

id H w= 0, (19) 

id tl w= -d^ w-2\w\ 2 w-w 2 w*. (20) 

The higher-order terms satisfy fist order ordinary differential equations with 
respect to t\. The spatial variable x\ is a parameter in these equations. The 
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solutions of these equation are uniquely denned by terms of the order of 1 in 
asymptotics as t\ — ► — oo. The asymptotics as t\ — > — oo is obtained by matching 

2 1 

w= d t2 u (xi,t 2 ) \t 2 =oh + 
w= -iL(2i\f\ 2 u)\ t2=0 + o(l), 



w = 7; d t 2 u (xi,t 2 )\ t2 =otl + u{ Xl ,0) + o(l). 



3 1 a9 1 / , \ i ,9.3, 

To determine the behaviour of the solution after resonance we need to calcu- 
late the asymptotics as t\ — > +oo of the coefficients for representation l|15fl. 
Calculations give 



w (xi,h) =u (xi,0) - if{x x ) 



ic x + eXP ff /2) + 0{tf) 



where c\ = (1 — i)\pH. 
Denote by 



2 

The function it) (xi, ti) has the asymptotics of the form 

2 . , 2 . . 2 . . , . exp(it?/2) _ 4 
to =tiw i(xi)+ io (a:i) + 5i(xi) ^ h 0{t 1 ), 

where 

2 ^2 1 i 1 i2 1 

w i = -»x! w o +| w | w; 

w (xi) = lirn^ ( J 1 w(x 1 ,6) + \w(x 1 ,6)\ 2 w {xi,0)]dB- tX 

gi(xx) — k\d Xl f + fc2|/| 2 /, ki and k 2 are constants. 

to (xi,ti) = i 2 w 2 (xi) +o(t 2 ), ti-^oo, 

where 

3 / 2 1 o 2 X o 2 \ 

w 2(^1) = i (dx x w i + 2 I w o I 2 »i + w w 1 ) > 

The domain of validity for (|15|l is defined by the same way as was shown in 
previous section. We require the following relation is valid 

e 5 p(ti, xi, e) — o(e 2 ), t\ — ► 00. 

The determined above behaviour of coefficients of asymptotics 1)150 gi ye the 
domain of validity for i|15|) 

M-^e- 1 or |i|<e" 2 . 



8 



4 Scattered wave 

In this section we construct the asymptotic solution of equation (JTJ after passage 
through the resonance. The leading-order term of the solution satisfies NLSE 
and depends on Xi,t% as well as before resonance. But this leading-order term 
is determined by another solution of NLSE which contains generally speaking 
another number of solitons. 

We construct the asymptotic solution of the form 

*(x,*,e) = e v (xi,t 2 ) + e 2 v (xi,t 2 ) + 

2 4 

e 2 A (t 2 ,^i)exp(iS'/e 2 ) + e 4 (Ai (t 2 , x t ) exp(iS/s 2 ) + 

A-i (t 2 ,£i)expHS/e 2 )). (21) 
Substitute this representation into (JTJ: 

e 2 (-S' A -/) exp(iS/e 2 ) + e 3 {d t2 v +d 2 Xl v+\l\ 2 l) + 

2 2 12 1^2 

e 4 (d t2 v +3^ « _|_2| « | 2 u + « v * + 
(-5' ii +a t2 1 +5^ 1 +2\ v | 2 1) exp(iS/e 2 ) + 
(S" l-i + w 2 1 *)exp(-^/e 2 )) = e 5 r(i2,a;i,e). 
Here r(^2,^i,e) depends on coefficients of the asymptotics (|21H . This depen- 

2 4 4 

dence is easy calculated. The coefficients A, Ai and A-i have singularity on 
the resonance curve. To determine the domain of validity of (|21[) we need to 
derive the explicit formula for r 

r(t 2 ,x 1 ,e) = 0(l + e\A\ 3 + s\n\e\+e 7 (\Ai \ 3 + \ A-i | 3 )). 

Collect the terms of the same order of small parameter and the same expo- 
nents. It yields the equations for coefficients of representation i|21|l . 

d t2 v +d 2 Xl l+\v\ 2 v= 0; (22) 

d t2 I +0 2 X1 v +2| v | 2 v + v 2 I * = (23) 

S' A= f; 
-S' Ai= -d t2 A -d 2 Xl A -2\ I | 2 A] 

s> L=-i 2 l*. 

Initial conditions for differential equations for v are obtained by matching. 
These conditions are evaluated on the resonance curve to = 0. 

v\ t2 = o =u(x 1 ,0) + (l-i)^f(xi); (24) 
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v | i2=0 =w (xi). (25) 

The residue part e 5 r(t2,xi,e) — o(e 2 ) as t 2 3> e. This condition is deter- 
mined the domain of validity for i|21|) . 

Formula l|24|) is connection formula for the leading-order term of the asymp- 
totic solution before and after the resonance. Additional term (1 — i)y/irf(xi) 
leads to changing of the solution after passage through the resonance. 



5 Numerical justification of asymptotic analysis 

In this section we justify our asymptotic formula Let us consider the pure 
soliton initial condition for equation 

_ 2V2ej 1 exp{-i2csx- 4(c 2 - V 2 )t e 2 } 
* (M ' £)|t = t0 " cosh^ex + 8cr 1 e 2 t + s) W 

According of our analytical results this initial condition leads to one soliton 
solution as the leading-order term of the asymptotic solution: 

l . . 2v / 2?7exp{-i2cxi - 4i(c 2 - r] 2 )t 2 } 

U (X\,t2) = -, r ■ 

cosh(2?7a;i + 8cr)t 2 + s) 

This soliton propagates up to the resonance curve t = 0. 

To annihilate this soliton on the resonance curve one may choose the specific 
form of the amplitude of the perturbation such that the left hand side of relation 
|@J equals zero: 

Q=i (<Bi,0) + (l-iWSF/(a:i). 

Hence _ 

-(1 + i) l -(l+i)y/2r]exp{-i2cx 1 } 

= o r u ^i' = 7= v75 ITT - 

2-^/Tr y 7r cosh(2?7a;i + s) 

To illustrate this by numerical simulations we choose e = 0.1, = 1, s = 0, c = 
0, t° = —200. Then the original equation 0) has the form 

idtV + a 2 * + |*| 2 * = 0.01/ exp{i0.005f 2 }. 

Initial condition is _ 

Tl 0.2V2 
* t=— 200 



cosh(0.2a;) 
and amplitude of the perturbation is 

-(1 + i) V2 



f 



\J~ii cosh(0.2a;) 



The numerical simulations of annihilation process for soliton of NLSE are pre- 
sented on the following figure. This justifies the formulas obtained above by 
matching method. 
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6 Summary 



We have studied the scattering process on the resonance and obtained the way 
for the control the parameters of scattered solitons. Using formula Q one 
can choose the form of the perturbation to obtain the scattered pattern of the 
given form. The process of the scattering on the local resonance is an universal 
phenomenon for waves propagate in dispersion media. We believe that our 
work can be useful for further understanding and description of solitary waves 
evolution under perturbation. 
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Figure 1: This picture shows the scattering process of one soliton to 
two solitons for equation l(T|. where amplitude of external force e = 0.1, 
S = 0.005i 2 , / = 2V2cosh" 1 (0.2x)+ 2^2 exp(0.2ia;) cosh _1 (0.2x - 5) + 
2v / 2exp(— 0.2ix) cosh _1 (0.2x + 5), initial data is pure soliton of NLSE at 
^|t=— 200 = 0.2^2 cosh" 1 (0.2a; ) , resonance curve is t = 0. 
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Figure 2: Annihilation of soliton. 
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